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ABSTRACT 

We investigate the form of the momentum distribution function for pro- 
tons and electrons in an advection-dominated accretion flow (ADAF). We 
show that for all accretion rates, Coulomb collisions are too inefficient to 
thermalize the protons. The proton distribution function is therefore deter- 
mined by the viscous heating mechanism, which is unknown. The electrons, 
however, can exchange energy quite efficiently through Coulomb collisions 
and the emission and absorption of synchrotron photons. We find that for 
accretion rates greater than ~ 10~ 3 of the Eddington accretion rate, the 
electrons have a thermal distribution throughout the accretion flow. 

For lower accretion rates, the electron distribution function is determined 
by the electron's source of heating, which is primarily adiabatic compression. 
Using the principle of adiabatic invariance, we show that an adiabatically 
compressed collisionless gas maintains a thermal distribution until the parti- 
cle energies become relativistic. We derive a new, non-thermal, distribution 
function which arises for relativistic energies and provide analytic formulae 
for the synchrotron radiation from this distribution. Finally, we discuss its 
implications for the emission spectra from ADAFs. 

Subject headings: accretion, accretion disks — black hole physics — radiation 
mechanisms: thermal, non-thermal — galaxies: elliptical and lenticular, cD 
— radio continuum: galaxies — Galaxy: general 



1. Introduction 

An advection-dominated accretion flow is a hot, optically thin, accretion flow with 
low radiative efficiency (Abramowicz et al. 1988; Narayan & Yi 1994, 1995a, 1995b; 
Abramowicz et al. 1995). Unlike standard thin disks (see Frank et al. 1992), where all 
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of the viscously generated energy is thermalized and radiated locally, ADAFs store most 
of the viscously generated energy as internal energy of the ions and advect it onto the 
central object. The gas in an ADAF is a two temperature plasma, with the ions being 
significantly hotter than the electrons (Shapiro et al. 1976; Rees et al. 1982). Since 
most of the viscously generated energy is assumed to heat the ions, and only a small 
fraction of this energy is transferred to the electrons via Coulomb collisions, the total 
energy radiated is much less than the total energy generated by viscosity (Rees et al. 
1982). 

The emission spectrum from an ADAF is determined by the cooling mechanisms 
of the protons and electrons. The protons cool very inefficiently, primarily through 
proton-proton collisions which create neutral pions. The pions decay into 7-rays with 
mean frequencies ~ 10 24 Hz (Mahadevan et al. 1997). The electrons, on the other 
hand, cool very efficiently through synchrotron, bremsstrahlung, and Compton processes. 
Detailed calculations of these cooling mechanisms lead to ADAF models which have been 
successfully applied to explain the observed radio to gamma-ray spectrum from a number 
of accreting black hole systems. These include solar mass black holes in X-ray binaries, 
e.g., A0620 and V404Cyg (Narayan, McClintock, & Yi 1996), as well as supermassive 
black holes at the centers of galaxies, e.g., Sgr A* (Narayan et al. 1995; Mahadevan et 
al. 1997) and NGC 4258 (Lasota et al. 1996). 

The shape and magnitude of an ADAF spectrum depends primarily on the mass 
accretion rate and, to a lesser extent, on the mass of the accreting object (Narayan 
& Yi 1995b; Narayan 1996). In determining the spectrum, however, all previous work 
has assumed that the electrons are thermal throughout the ADAF. This is a crucial 
assumption, since the emergent spectrum would be considerably different for different 
electron distribution functions. In the present paper, we investigate the validity of this 
assumption. 

We treat thermalization within the ADAF paradigm. We assume that the only phys- 
ical interactions among particles are Coulomb collisions and self-absorbed synchrotron 
radiation. We do not consider the possible non-thermal heating of electrons by wave- 
particle interactions, magnetic reconnection, or collective plasma effects. The interac- 
tions considered here therefore incorporate only the well understood physical processes 
in ADAFs. 

The outline of the paper is as follows. In the following section (§2) we present 
the necessary background equations and consider proton thermalization by Coulomb 
collisions. In §3 we investigate the form of the electron distribution function, in particular 
for suprathermal electron energies, and determine the mass accretion rates for which 
electrons are thermalized by Coulomb collisions or self-absorbed synchrotron radiation. 
In §4 we consider the form of the electron distribution function when thermalization is 
inefficient, which depends on the interplay between the heating and cooling mechanisms. 
We give a new (non-thermal) distribution function, valid when the electrons compress 
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adiabatically to relativistic energies. In §5 we give the synchrotron spectrum from this 
new distribution function and in §6 we discuss the implications and possible applications 
of our results. 



2. General Relations 

2.1. Self-similar Flow Equations 

The parameters which determine the structure of an ADAF are the viscosity pa- 
rameter, a (Shakura & Sunyaev 1973), the ratio of the gas pressure to the total pressure, 
/3 a dv, the mass of the central object, M = m M , where m is the mass in solar mass 
units, and the accretion rate, M = rh M^d, where m is the accretion rate in Eddington 
units (MEdd = 1-39 x 10 18 m g s _1 ). Typical values for a and /3 a d v which have been suc- 
cessfully applied to observed systems are 0.3 and 0.5, respectively (Narayan et al. 1996, 
1997). /3 a dv = 0.5 corresponds to equipartition between magnetic and gas pressure, as 
suggested by the work of Balbus and Hawley (1991). 

Unlike standard thin accretion disks, an ADAF is well approximated by a series 
of concentric spherical shells with the properties of the gas varying as a function of 
radius (Narayan & Yi 1995a). The self similar solution of Narayan and Yi (1994, 1995a) 
provides reasonably accurate analytical estimates of the properties of the accretion flow. 
For the present discussion, the quantities of interest are 

c 

P {r) - 0.09 r-\ 
n e {r) ~ 6.3 x 10 19 a _1 m _1 mr 3 / 2 cm -3 , 

B = 7.7 x 10 8 a' 1 ' 2 (^f^) V2 m-^fn 1 ' 2 r" 5 / 4 Gauss, (1) 

where v(r)/c is the radial velocity in units of the velocity of light, 9 p (r) = kT p /m p c 2 is 
the dimensionless proton temperature, n e (r) is the number density of electrons, B(r) is 
the magnetic field strength and r = R/Rs is the radius in Schwarzschild units (Rs = 
2.95 x 10 5 m cm). In equation flU), the fraction, /, of the viscously dissipated energy that 
is carried inward by the accreting gas is taken to be ~ 1 (Narayan et al. 1996). 
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2.2. Timescales 



For the electrons or protons in an ADAF to remain thermal as they accrete, they 
must have sufficient time to redistribute their kinetic energy by various processes. We 
therefore take thermalization by Coulomb collisions to be efficient when the timescale 
for collisions between the electrons, t ee , or between the protons, t pp , is shorter than the 
accretion time, t a . Since the protons are non-relativistic (8 P < 0.1) their thermalization 
time is (Spitzer, 1962) 

rip(JTC In A \m e J ' 

~ 9.2 x 10~ 3 a (^^j ' mm- 1 s, (2) 

where n p = n e is the proton number density, <jt is the Thomson cross-section, In A ph 20 
is the Coulomb logarithm, and we have used equation (|l|) for the proton temperature 
profile. The time required for the protons and electrons to reach thermal equilibrium is 
(Spitzer 1962) 

(27T) 1 / 2 (m p \ /2 
tep ~ 2n e cr T c\nA {^~J ^ + 9p) ' 

~ 9.3 x 10- 5 a9 3 J 2 mm' 1 r 3/2 s. (3) 

So long as the electrons are non-relativistic their timescale for thermalization is (Spitzer 
1962) 



(2vr) 1/2 



Q 3/2 



n e aTC In A 

~ 1.0 x 10" 7 a mm' 1 r 3/2 6 3 J 2 s. (4) 

For thermalization to occur before the gas falls into the central black hole, the above 
timescales must be less than the accretion time of the gas, which is 

t a = /4^y - 1-8 x 10" 5 a" 1 mr 3/2 s. (5) 

From these timescales we can readily estimate the critical accretion rate, m crit , 
above which the ADAF solution ceases to exist. For sufficiently large mass accretion 
rates, the electrons and protons become thermally well coupled and so the assumption 
of a two temperature plasma is no longer valid; the two temperature ADAF therefore no 
longer exists. The accretion rate above which this occurs is obtained by setting t cp ~ t a , 
which gives 

m crit ~ 0.45 a 2 , (6) 

where we have set 9 e ~ 0.2, which is valid for high m systems (Mahadevan 1997). This is 
comparable to the critical accretion rates obtained by independent arguments (Narayan 
& Yi 1995b; see also Mahadevan 1997). 
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2.3. Proton Thermalization 



To determine the accretion rates at which the protons have sufficient time to ther- 
malize by Coulomb collisions, we set t pp ^ t a . This gives the requirement 

m> 500a 2 ^^'V 3 ' 2 . (7) 

Since the maximal accretion rate for an ADAF is m cri t, the radius above which the 
protons can thermalize is 

For a given rh, the protons are unable to redistribute their energy through Coulomb 
collisions for r < r t h and so the proton distribution function at these radii is primarily 
determined by the viscous heating mechanism. Observational probes of the proton 
distribution function would therefore provide valuable information about the details 
of the viscous heating mechanism. Mahadevan et al. (1997) have shown that the 7- 
ray spectrum from pions created by proton collisions in ADAFs can be used as a direct 
probe of the proton distribution function, and can therefore provide a means of studying 
viscosity in hot accretion flows. 



3. The Electron Distribution Function in ADAFs 



In this section we analyze the form of the electron distribution function in ADAFs. 
We show that for moderately high accretion rates, Coulomb collisions and synchrotron 
self-absorption can efficiently thermalize the electrons, even in the high energy tail of 
the distribution function. At lower accretion rates, however, these processes are too 
inefficient, and the form of the distribution function depends on the dominant heating 
and cooling mechanisms of the electrons. This is discussed in §4. 



3.1. Coulomb Collisions 



The accretion rate, m t h, above which the electrons can thermalize via Coulomb 
collisions is obtained by setting the timescale for electron-electron collisions, t ec , equal 
to the accretion time, t a , which gives 

m th « 3.8 x l(r 4 T 9 3/: V, (9) 
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where T 9 is the electron temperature in units of 10 9 K. For rh <: m t h, electrons with 
energies ~ kT e are thermal. Since most of the synchrotron emission from ADAFs comes 
from the high energy tail of the distribution (§4), it is also necessary to determine whether 
electrons with energies 3> kT e have enough time to thermalize. To investigate this, we 
use the Fokker-Planck equation which, for an isotropic energy distribution function, 
takes the form 



dt 



d_ 

$7 




n(7,t) 



+ 



1 d 2 



2 



dt 



(10) 



where n(j,t) is the distribution function, dj/dt represents systematic acceleration, and 
d{/\^) 2 /dt represents stochastic acceleration. The systematic acceleration is negative 
for particles with 7 ^ 9 e , and positive for 7 ^ 9 e , which corresponds to slowing down 
(speeding up) particles with energies large (small) compared to the mean energy of the 
plasma. The stochastic acceleration is always positive and corresponds to diffusion in 
energy space. 

The timescale for an electron with a Lorentz factor 7 to thermalize with a back- 
ground plasma of temperature 9 e is given roughly by the timescale on which the electron's 
energy e-folds. This ensures that the electron interacts significantly with the plasma. 
The resulting timescale depends on 7 and 9 e , and is obtained by determining whether 
systematic or stochastic acceleration dominates the energy changes of the electron. 

The systematic acceleration of an electron with Lorentz factor 7, interacting with 
a thermal plasma of temperature 9 e , is given approximately by (Dermer & Liang 1989) 



g?7 
~dt 



— - n e OTC In A 



1 



A(9 e 



(11) 



where A(9 e ) = K 2 {\-/9 e ) / Ki{l/ 9 e ) and K n is the modified Bessel function of order n. 
The timescale on which the electron's energy e-folds is therefore 

(12) 



t e 



7 



\d\wy/dt\ 3 n e o"TclnA 



1 



A(9 e 



Equations (|11]) and ( |T2"D are good analytical approximations to detailed numerical cal- 
culations provided that 9 e <: 0.3 and 7 2 (Dermer & Liang 1989). 

The stochastic acceleration is determined by substituting equation ([TT]) into the 
Fokker-Planck equation, setting dn/dt = (as is required for a thermal plasma) and 
solving the second order differential equation. The two integration constants in the 
resulting solution are determined by (1) requiring no net flux of electrons in energy 
space (since there are no sinks or sources of electrons) and (2) imposing the physical 
requirement that the stochastic acceleration not increase exponentially for large 7. This 
gives 



dt 



3n e <7Tcln A 



9 r , 



A(9 e ) 7 



e 

+ ~Z6 e + 



01 



->ze e 



(13) 
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where z$ e = 2[6 e / A(9 e )] — 1. The resulting diffusion time is 



6 diff 



/ 1 d(A 7 ) 2 

1 « 2 dt 



r 



el 



A(6 e 



H ^ e + -o^tfe 

7 7 



(14) 



7^ dt 3 n e GTC In A 

Since equations ([X3| ) and ( |I4"D are derived from equation (|Tl"|), they are also only valid 
for 6 e <: 0.3 and 7^2. We note that the stochastic acceleration derived here does not 
rise exponentially as indicated in Dermer & Liang (1989; cf. their eq. [14] and figure 2). 
More recent detailed calculations by Nayakshin & Melia (1997) agree with the analytic 
approximation for the stochastic acceleration given here. 



Setting (t sys ,tdifr) < and using equations (H) and 
required for thermalization, 

< 1.1 x ur :! 



gives the accretion rates 



Q'.~ 



sys 
c 



m \ 
a 2 J diff 



> 



7 



7.3 x 1(T 4 7 2 



A(e e ] 



A(6 e 



+ —Z6e + 

7 



(15) 



For electrons with energies comparable to the mean energy of the plasma, 7 ~ 9 e and 
equation ([12]) shows that tl ys — > 00; the electron, on average, is not heated by direct 
acceleration, but instead diffuses in energy space (as determined by the stochastic term). 
For large 7, tg ys < t^ iS , and so electrons with energies much larger than the mean energy 
of the plasma, 7>0 C , first lose energy by direct cooling and then relax by diffusion. In 
order to assess the efficiency of thermalization, we therefore take the minimum of the 
systematic and diffusion timescales as the relevant timescale against which to compare 
the accretion time. 

For a plasma at temperature 8 e , an electron with Lorentz factor 7 can thermalize via 
Coulomb collisions if the accretion rate is greater than the smaller of the two accretion 
rates given in equation (|i~5"l). The dashed lines in figure 1 show this accretion rate as a 
function of the Lorentz factor 7, for two values of 9 e . For accretion rates below these 
curves, the electrons do not have time to interact with one another via collisions and 
are "frozen" into their initial distribution. Unless other thermalization processes are im- 
portant, this implies that the form of the distribution function can deviate substantially 
from a Maxwellian. We show below that because the synchrotron photons in ADAFs are 
highly self-absorbed, the electrons are able to thermalize at lower accretion rates (and 
higher 7) by interacting via the exchange of self-absorbed synchrotron photons. 



3.2. Electron Thermalization Through Synchrotron Self- Absorption 

Synchrotron self-absorption has been treated in the context of finding steady state 
power-law electron distributions (McCray 1969), heating low energy electrons through 
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self-absorption (Ghisellini et al. 1988), and as a thermalizing mechanism for low energy 
electrons (Ghisellini & Svensson 1990). In particular, Ghisellini and Svensson (1990) 
have shown that self-absorption of synchrotron photons leads to a thermal distribution 
of electrons over a large range of 7 in a few synchrotron cooling times. The basic ther- 
malization mechanism is energy exchange via the absorption and emission of synchrotron 
photons, which can be important even when the plasma is effectively collisionless. 

In keeping with the analysis of the previous section we discuss thermalization in- 
duced by synchrotron self-absorption using the Fokker-Planck equation. Instead of 
attempting to solve the complete integrodifferential equation (eq. |IJ), we check for con- 



sistency by assuming that the (thermal) synchrotron radiation is highly self-absorbed, 
and then determining whether the thermalization timescales obtained from the Fokker- 
Planck coefficients are faster than the accretion time. 

The kinetic equation governing the electron distribution function can be written 
as a Fokker-Planck equation with coefficients given by (McCray 1969; Ghisellini et al. 
1988) 

d"i 1 r°° Id 



dt m e c 2 Jo ' 7P (?7 

d(A 7 ) 



dt 



2C(7). (17) 



1 z 100 I( u) 
C(7) = 7r-rW -y-M-r)^. 



where 

nM = 

2m 2 c 2 Jo v 2 

7) is the cyclo-synchrotron emissivity (ergs s" 1 Hz -1 ), p = 7/3 = (7 2 — l) 1 / 2 is 
the electron's momentum in units of m e c, and I(y) is the intensity of the background 
photons (ergs cm -2 s _1 Hz -1 ). 

To evaluate C ('-/), we use three properties of synchrotron radiation in ADAFs. First, 
the electrons which are responsible for most of the synchrotron emission are from the 
high energy tail of the distribution. By multiplying the relativistic Maxwell-Boltzmann 
distribution by the synchrotron emissivity and using the method of steepest descent, 
we can determine the Lorentz factor (7 max ) of the electron responsible for most of the 
emission. This gives 

7max = 2 1 /%x\ / I 3 ~ lO0 e , (19) 

where %~ 500 and is a very weak function of rh (Mahadevan 1997) . Large values of xm , 
as is the case in ADAFs, imply that most of the observed synchrotron radiation originates 
in the exponential tail of the single particle emissivity. Since the electrons responsible 
for most of the synchrotron emission are relativistic (eq. we can evaluate C(7) by 

using the synchrotron formula for j(u,j) (Rybicki & Lightman 1979) 

j(y,j) = S sina p F ( — 
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3 1 / 2 (2 7 r)e 2 z/ B 3 2 . e_B 

S = , Vc = 7:7 v B sin Op, z/ B = , (20) 

c 2 Z7rm e c 

where a p is the angle the electron makes with the magnetic field, F(x) is an integral over 
modified Bessel functions (Rybicki & Lightman 1979, eq. [6.31c]), and vb = 2.8 x 10 6 B 
Hz is the cyclotron frequency. 

Second, an electron with Lorentz factor 7 max emits most of its radiation at its 
critical frequency u c ~ 1.5 J^J^b- This radiation is, however, highly self-absorbed and 
becomes optically thin only at a frequency v t = 1.58 2 xmVb ~ 4.7z/ c (Narayan & Yi 
1995b; Mahadevan 1997). We can therefore set u t ^> v c to good approximation. 

Finally, since the plasma is highly self-absorbed, the radiation field, liv)-, is Raleigh- 
Jeans at a temperature 9 e up to the frequency u t , where the plasma becomes optically 
thin (see Narayan & Yi 1995b; Mahadevan 1997). This gives 

(2v*m e e e 0<u<u t , 
1 ; I v > v t . K ' 



Combining these properties, equation (|l8j) becomes 

3 1/2 e 3 B r u t/ v c 

= n~ / F(x)dx, 

m^r jo 

~ 1.2 x l(T 9 £ 2 e 7 2 , 

= C.BX^ 2 , (22) 

where C\ is a constant as defined above and we have averaged over all pitch angles. 

Equations (1TB) and (O) are now easily evaluated to give the systematic and diffusion 
timescales, 



7 



sys \d\n-f/dt\ dB^ 

f 1 / 1 d ^ - 1 (23) 

Note that in the cyclotron limit, the systematic and stochastic coefficients can be eval- 
uated by recognizing that j(v,y) oc 5{y — Ub), which yields 

dt 

d -^- ~ dtfS,^- (24) 
at 

In the limit /3 — > 0, the stochastic term is zero, and the electrons only heat up by 
absorbing radiation from the background field (Ghisellini & Svensson 1990). Since we 
are concerned with thermalization of high energy electrons, we neglect the cyclotron 
limit in our discussion. 
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Setting (ts ys ,tdiff) < ^ai an d using equations (]!]) and (|5|), gives the accretion rates at 
which thermalization by self-absorption of synchrotron photons is efficient 

^ 7.6 x 10~ 5 - 

OT J sys 7 

a 3.8x10-1 fm*) " r. (25) 

For a plasma of temperature 9 e , an electron with Lorentz factor 7 can thermalize through 
synchrotron self-absorption if the accretion rate is greater than the smaller of the two 
accretion rates given in equation (p5|) . The solid lines in figure 1 show this accretion rate 
as a function of the Lorentz factor 7 for two values of 9 e at r = 10. For fixed 9 e , these 




curves can be scaled linearly with r for other radii (cf. eq. p5|). As Figure 1 indicates, 
synchrotron self-absorption is significantly more efficient than Coulomb collisions at 
thermalizing high energy electrons. 

Note that unlike electron thermalization, thermalization by synchrotron self-absorption 
depends explicitly on the radius. The magnetic field is weaker at larger radii which, for 
a fixed temperature, decreases the amount of synchrotron radiation. This decreases the 
amount of self-absorption in the plasma and so thermalization is less efficient, as is clear 
from equation (p5|). 



The systematic term in equation fl2"o] ) is the dominant thermalizing mechanism for 
electrons with 7 ~ 7max- This term can be rewritten to give a condition on the radius 
below which the plasma is thermal, 

r t W.3xl0 4 7 g), (26) 

where we have taken 7 3> 9 e . For a given 7 and m/a 2 , the electrons are thermal for all 
radii r ^ r^ h . To determine the accretion rates above which electrons with 7 ~ 7max are 
thermal throughout the flow, we set 7 ~ 106^ and rf h ~ 1000 in equation (|26|), which 
gives 

rh 

In the last equality, we have set 9 e ~ 0.2, as is valid for high accretion rates (Mahade- 
van 1997). To determine the accretion rates at which the electrons cannot thermalize 
anywhere in the ADAF, we set r£ h ~ 10 and obtain 

Tfl 

— <> ~ lO" 4 , (28) 

cr 

where we have set 9 e ~ 1, which is valid at small radii for low accretion rates (§4.2). For 
these low accretion rates, the electrons responsible for most of the synchrotron emission 
are unable to thermalize throughout the accretion flow, and the assumption of thermal 
synchrotron radiation is no longer valid. The form of the electron distribution function 



— > 8 x lO^J 1 ~ 4 x 10~ 2 , (27) 
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and the resulting synchrotron spectrum at these low accretion rates are determined in 
§4 and §5. 

For accretion rates between these extremes, 10 -4 ^ raj a 2 ^4x 10~ 2 , the electrons 
will be thermal at small radii where synchrotron self-absorption is efficient, but will not 
necessarily be thermal at larger radii. The implications of this are considered in §6. 

To conclude this section, we note that the above analysis breaks down for electrons 
with large Lorentz factors (7 ^> 7 ma x)- This is because the radiation produced from 
these electrons is not self-absorbed, and escapes the plasma freely. The critical Lorentz 
factor, 7 cr it, at which this occurs is estimated by setting the electron's critical frequency 
equal to the frequency at which the plasma becomes optically thin, v t = l.hj^yB- This 
gives 

7crit = 6Um 2 , (29) 

which is always greater than 7 max for ADAFs since 7 cr i t = 2~ 1 / 3 x] v ( 6 7 max ^ 27 max .[] For 
electrons with 7 > 7 cr ; t the stochastic Fokker-Planck coefficient (cf. eq.[]17| 



is zero, 



since there is no absorption, and the systematic term is due solely to optically thin 
synchrotron cooling (cf. eq.[]T6 with C{^f) = 0). The electron distribution for 7 > 7 crit 



is therefore not necessarily thermal. This, however, has little effect on the synchrotron 
spectrum since most of the emission comes from electrons with 7 « 7 max < 7 cri f 



4. Heating and Cooling of Electrons 



When thermalization by Coulomb collisions and synchrotron self-absorption is inef- 
ficient, the evolution of the electron distribution function is determined by the interplay 
between the dominant heating and cooling mechanisms. These can be determined by 
considering the electron energy equation 

ds de , _ 

P ~dR = nV ~dR ~ q = q q ' ( 30 ) 

q ~ kIV dR' 
q e+ = q ie + q v , 

where s is the entropy of the electrons per unit mass of the gas, e is the internal energy 
of the electrons per unit mass, q c is the compressive heating (or cooling) rate per unit 
volume, and q~ is the energy loss due to radiative cooling. The total external heating 



however for accretion rates below ~ 10 9 the plasma is no longer self-absorbed and there is no 
thermalization by synchrotron self-absorption (7 cr it < 7max since xm < !)■ 
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of the gas, q c+ , is a sum of the heating via Coulomb collisions with the hotter protons, 
q ie , and direct viscous heating, q v . 

Much of the previous work on ADAFs has neglected the compressive heating of 
electrons and has assumed that the electron temperature profile is given by solving 
q- = q c+ (e.g., Narayan & Yi, 1995b). We show below that this is likely to be a 
poor approximation except for systems near the critical accretion rate. This point has 
already been made by Nakamura et al. (1997), who first recognized the importance of 
compressive heating in determining the electron temperature profile. 

As long as q c+ <C q c and q~ <C q c , the right hand side of equation (|30|) is small 
compared to each term on the left hand side. All of the compressive heating therefore 
goes into changing the internal energy of the gas, with nearly constant entropy. We 
now show that both of these conditions are satisfied at low accretion rates, so that the 
electrons in ADAFs compress adiabatically as they accrete onto the central object. 



4.1. Heating 

The dominant heating mechanism for the electrons is obtained by comparing the 
relative magnitudes of q c , q m , and q v . Using equations ([!]) and ([3(]), the compressive 
heating rate is 

g c ~5x 10 17 mm~ 2 T 9 r~ 3 erg cm" 3 s _1 . (31) 

The energy transfer rate between protons and electrons via Coulomb collisions is (Step- 
ney & Gilbert 1983) 

q ie ~ 2.05 x 10- 17 n e n 4 (T, - T e )T^ 2 

= 2.5 x 10 21 a- 2 m 2 m- 2 T 9 3/2 r- 4 erg cm" 3 s" 1 , (32) 

where we have assumed that T e < 10 10 K so that a non-relativistic formula is reasonably 
accurate. The last equality follows from using equation ([!]), and setting T e «C T i7 n e = rij. 

The viscous heating of the electrons is determined by allowing for a fraction 5 ~ 
m e /m p ~ 1/2000 of the viscously generated energy to be transferred directly to the 
electrons, so that (Mahadevan 1997) 

q v ~ 3 x 10 17 mm _2 r~ 4 <5 o erg cm -3 s -1 , (33) 

where 8 = 20005. Both compressive and viscous heating are proportional to the accre- 
tion rate and their ratio is given by 



q e /q v ~2 6; 1 T 9 r, 



(34) 
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which is independent of the global parameters of the ADAF. In the inner regions of the 
accretion flow (r ~ 3), T g ~ 5, while for larger radii the temperature decreases roughly 
as r -2 / 3 (§4.2). Compressive heating of the electrons is therefore always more important 
than direct viscous heating provided that 8 is not substantially larger than 10 -3 . 

Coulomb heating of the electrons by the protons is, however, a two body process 
and becomes less important for small m. From the expressions for q c and q ie given above, 
it follows that compressive heating dominates over Coulomb heating provided that 

m < 2 x l(TVT 9 5/2 r. (35) 



Comparing this with equations (|T5|) and (|2~5 ) shows that for low m, where both Coulomb 



collisions and synchrotron self-absorption are unable to thermalize the plasma, compres- 
sive heating is the dominant heating mechanism. 



4.2. Cooling 



Synchrotron cooling can be an important energy loss mechanism in ADAFs because 
of the assumption of near equipartition magnetic fields. In fact, for the low mass ac- 
cretion rates of interest here, synchrotron radiation is the dominant source of cooling 
(Narayan Sz Yi, 1995b; Mahadevan 1997). In §2 we showed that the thermalization 
time through synchrotron self-absorption is of order the electron cooling time and that 
self-absorption is inefficient at thermalizing the electrons for m ^ 10 a 2 (cf. Figure 1). 
For these low accretion rates electrons therefore cannot radiate a significant fraction of 
their energy in an accretion time. Since the dominant source of the electron's energy is 
compressive heating (q c ) this implies that q~ <C q c . 

This argument relies on the accuracy of our estimate of the electron temperature in 
the accretion flow, since higher temperatures would lead to more efficient synchrotron 
cooling. Since previous work has neglected compressive heating, however, we must verify 
that electron temperatures of ~ 10 9 K are still valid when the electron temperature 
profile is determined by adiabatic compression. 

For adiabatic compression, the temperature profile is given by T oc p r_1 , where V is 
the ratio of the gas specific heats and p oc r~ 3 / 2 . Since the internal energy of the electron 
gas includes a contribution from the magnetic field energy, T = (8 — 3/3 a d v )/(6 — 3/9 a d v ) 
for non-relativistic electrons (Esin 1997). For an equipartition magnetic field, V = 13/9, 
which gives a temperature profile of T oc r~ 2 / 3 . T = 4/3 is still valid in the relativistic 
limit, in which case T oc r~ a5 . For T e »s T p at r m 10 4 , the above scalings imply that 
the electron temperature near the last stable orbit is ~ 5 x 10 9 K, which is comparable 
with previous estimates. 
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4.3. The Distribution Function of an Adiabatically Compressed Gas 

For m ^ 10 _4 a 2 , Coulomb collisions and synchrotron self-absorption are unable 
to thermalize the electrons. The electrons therefore do not interact with one another 
and can be well approximated as a collisionless gas. Furthermore, at these low accretion 
rates, the electrons compress nearly adiabatically as they accrete onto the central object. 

An adiabatic change is one in which the properties of a system vary slowly compared 
to the characteristic timescales of the system. During such a change, the system's phase 
space volume is a conserved quantity (Landau & Lifshitz, Mechanics §49). We use this 
property to determine the evolution of the distribution function of a collisionless electron 
gas subjected to adiabatic changes in its volume. We emphasize that the discussion 
which follows is not specific to an ADAF, but applies to any collisionless plasma which 
undergoes adiabatic changes. 

The adiabatic compression of a collisionless gas can be viewed as the independent 
adiabatic compression of each particle in the gas, since there is no coupling between the 
particles. By adiabatic invariance, the quantity p 3 p~ x is conserved for each particle in 
a monatomic gas, where p is the momentum of the particle and p is the density of the 
gas.0 The evolution of the distribution function is therefore given by the transformation 
p — > ap, where a = (pf/pi) 1 / 3 , and p^f\ denotes the initial (final) gas density. We 
emphasize that this transformation is valid for all momenta p, both non-relativistic and 
relativistic. 

To determine how the distribution function of a non-relativistic thermal gas of tem- 
perature 8 e changes under an adiabatic compression, consider the Maxwell-Boltzmann 
distribution function 

n{p)dp = (2 Jw/2 exp {-^he) dp - (36) 

For now, assume that the compression factor, a, is such that the particle momenta after 
compression remain non-relativistic. Substitution of the transformation p — > ap into 
equation (|36| ) shows that the form of the distribution function remains the same. The 
gas is still thermal, but is now characterized by a temperature 6' e = a 2 9 e = (pf / Pi) 2 ^ 3 e , 
which is the usual result for the adiabatic compression of a F = 5/3 gas. A thermal 
distribution is therefore maintained even though the gas is collisionless. 

Similarly, for a relativistic thermal gas characterized by the Maxwell-Boltzmann 
distribution function, 

n{p)dp = ^mktm exp (~^) dp (37) 



2 p 3 p 1 is also the adiabatic invariant for particles moving in an isotropically tangled, flux frozen, 
magnetic field. 
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setting p —> ap and assuming that the final momenta are relativistic shows that the gas 
remains thermal with a final temperature 6' e = a9 e = (pf / ' pi) l ^6 e . This is the usual 
result for the adiabatic compression of a V = 4/3 gas. 

Under adiabatic compression, an initially non-relativistic (relativistic) thermal gas 
therefore remains a non-relativistic (relativistic) thermal gas, as long as the post com- 
pression momenta remain non-relativistic (relativistic). Consider, however, a non- 
relativistic thermal gas (eq. |j3~6|)) which is compressed by a sufficient amount that 



the final momenta of the particles are now relativistic. Since the transformation p — > ap 
is valid regardless of the magnitude of the momenta, the post compression distribution 
function is given by 

Anp 2 ( p 2 \ 

Hip)dp = (2WcW e )3/ 2 ex P (- W ^ e ) dp - (38) 

This has the form of a non-relativistic Maxwellian, even though the particle energies 
can be highly relativistic. Therefore a collisionless gas which is adiabatically compressed 
across the non-relativistic/relativistic boundary does not maintain a thermal distribu- 
tion.Q In particular, the relativistic high energy tail of the adiabatically compressed 
distribution function falls off as a Gaussian rather than as an exponential (cf . eqs. [|38| . 

[0)- 

In order to assess the implications of this new distribution function, we compare 
it with a thermal distribution function which has the same average energy (the "best 
fit" thermal distribution). The average energy of an adiabatically compressed gas in the 
relativistic limit is (in units of m e c 2 ) 

4:71 

(7) = 



(2vr a 2 e e f' 2 Jo 
4 




e\' 2 . (39) 



(27^/2 

The average energy of the gas increases as a in the relativistic regime, rather than as a 2 
as in the non-relativistic regime. Since the average energy of a thermal relativistic gas 
is (7) = 39m b, the "best fit" thermal gas has 

6 M B=(^j 12 ae\' 2 . (40) 



Figure 2a shows the distribution function for a gas which is adiabatically com- 
pressed by a factor of a = 20, from an initially thermal distribution with 6 e = 0.01. 
Superimposed in the figure is the distribution function for a thermal gas which has the 



3 Note that for photons, which are massless, there is no such boundary, and so adiabatic changes 
always maintain a thermal distribution, a result which is, of course, well known in cosmology. 
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same mean energy as the adiabatically compressed gas (9mb — 1). As expected, the 
high energy tail of the adiabatically compressed gas is strongly suppressed with respect 
to the thermal gas. 

A relativistic collisionless gas that adiabatically expands across the relativistic/non- 
relativistic boundary will also not maintain a thermal distribution. Fig 2b shows the 
distribution function for a gas which adiabatically expands by a factor of a = 0.01 
from an initially thermal distribution with 9 e = 10. Superimposed in the figure is the 
distribution function for a thermal gas which has the same mean energy (9mb — 0.04). 
The high energy tail of the adiabatically expanded gas is overpopulated with respect to 
the thermal gas since the distribution function falls of as an exponential, rather than as 
a Gaussian, in the non-relativistic limit. 

Finally, we note that it is not necessarily valid to assume that the electrons have an 
adiabatically compressed distribution function just because thermalization is inefficient. 
This is because the timescale for thermalization as determined in §3 is not the same as 
the timescale on which an initially non-thermal distribution function changes substan- 
tially. As we discuss in Appendix C, it is only strictly true that the electrons have the 
adiabatically compressed distribution function out to 7 ~ 7 max if the accretion rate is 



~ 10 times smaller than the thermalization accretion rate given in §3 (cf. eqs. p5|). 
However, to accurately determine the form of the distribution function in this regime 
requires solving the complete Fokker-Planck equation, which is beyond the scope of 
this paper. In what follows, we assume that, for this narrow range in accretion rates, 
the electrons are well approximated as having an adiabatically compressed distribution 
function. 



5. The Synchrotron Spectrum from an Adiabatically Compressed Gas 

For an adiabatically compressed gas the large decrease in the number of particles 
at high energies has important consequences for the synchrotron emission from ADAFs. 
This is because most of the synchrotron emission originates from high energy electrons 
in the tail of the Maxwell-Boltzmann distribution. A lack of these high energy electrons 
will substantially reduce the synchrotron luminosity. 

To determine the synchrotron emission from this new distribution, we first calculate 
the synchrotron source function, S u , and then calculate the emergent intensity. In 
the present analysis, we assume that all of the electrons are relativistic so that the 
distribution function (eq. p8| ) takes the form 

4vr / 7 2 \ 

n(7) dl = (W^ 7 6XP [-2*% ) dr (41) 
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Using equation (|20|), the emissivity, 6(p)ac, an d the absorption coefficient, a(u) A c, are 
given by 

e{v) AG dv = J j(v, 7) n(j)dj, 

= E ^-I'(x' M ,2)diy ergs s" 1 Hz" 1 , (42) 



a{v) AG dv = -- — - — - I ' Jpj(v,l)S- 
aiim e v A J 07 



IP 



d7, 



Ei { J 9 ' )3/2 I'(x'm,3) dv cm" 1 , (43) 



where p — > m e ^c in the last equality. We have defined 
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and 



to maintain a similarity between the equations presented here and those for thermal 
synchrotron emission (Pacholczyk 1970; see also Mahadevan et al. 1996). I'(x' M ,n) is 
I(x' M ,n) integrated over all particle directions. The asymptotic expansions for I(x' M ,n) 
and I'(x' M ,n) at large and small x' M are given in Appendix A. 

The source function is determined using the limiting expressions for I'(x,n) (Ap- 
pendix A), which gives 

q e(u)AC j 2^my (a 2 ^) 1 / 2 x' M « 1, 

v 47ra(z/) AC \ 2 1 / 2 3 1 / 4 m e 4 / V/ 4 (a 2 # e ) 3 / 4 x' M > 1. 1 ' 

The source function for a thermal gas in the Raleigh- Jeans limit is S v = 2m e v 2 8 e . For 
low frequencies, x' M <C 1, the source function for the adiabatically compressed gas has 
the same frequency dependence as that of a thermal gas, but the normalization is larger 
by a factor of 3^/^T/A ~ 1.33 because there is an excess of low energy electrons (see 
Figure 2a). For large frequencies, x' M 3> 1, however, the frequency dependence of the 
source function for the adiabatically compressed gas is different than that of a thermal 
gas, and varies as z/ 7 / 4 . 

The emissivity for an adiabatically compressed gas at large frequencies is given by 

e(v) Ac dv = 1.6 x 1(T 28 u B x' M 3/4 exp (-2x' M 1/2 ) dv, ergs s" 1 Hz" 1 , (46) 

where we have used the expression for I'(x' M , n) for x' M 3> 1 (Appendix A). For a thermal 
gas, the corresponding emissivity is (Mahadevan et al. 1996, eqs. [11], [A10]), 

e(z/) Th dv = 3.4 x 1(T 28 x% 6 exp (-1.8899^ 3 ) du, ergs s" 1 Hz" 1 , (47) 
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where 

XM = W e ' 

The primary difference in the emissivity is the frequency dependence in the exponential. 
The exponential in the thermal spectrum falls off as z/ 1 / 3 , while it falls more rapidly, as 
z/ 1 / 2 , for an adiabatically compressed gas. This is expected because there is a lack of 
high energy particles in the adiabatically compressed gas as compared to a thermal gas 
with the same average energy. 

Figure 3a compares the total synchrotron luminosity from a sphere of radius R for 
a thermal and an adiabatically compressed electron gas with the same average energy 
(9 e = 0.5 for the thermal gas). We set R ~ 7 x 10 11 cm, B ~ 200 Gauss, and n e ~ 10 10 
cm -3 , which are the parameters at ten Schwarzschild radii for a 2 x 10 6 solar mass black 
hole accreting at rh ~ 10~ 5 . The solid line shows the spectrum for the adiabatically 
compressed gas, the dashed line shows the spectrum for the thermal gas, and the dotted 
line shows the optically thin emission that would result in each case were the plasma not 
self-absorbed. The adiabatically compressed gas has a significantly smaller synchrotron 
luminosity and becomes optically thin at a much lower frequency than the thermal gas. 
The adiabatically compressed gas is, however, still highly self-absorbed, thus validating 
our assumption that most of the emission comes from the exponential tail of the single 
particle emissivity. 

Since the electron number density and magnetic field strength vary with radius, 
the total synchrotron spectrum from an ADAF is obtained by summing the individual 
spectrum from each radius in the accretion flow (e.g., Narayan & Yi, 1995b). An example 
of this is shown in Figure 3b for an adiabatically compressed gas (solid line) and a thermal 
gas (dashed line), taking m = 2 x 10 6 , rh = 10~ 5 , and using equation (|l|) for the number 
density and magnetic field as a function of radius. In this calculation, we have assumed 
for simplicity, and for comparison with previous work, that the internal energy of the 
electrons is constant throughout the accretion flow (9 e = 0.5 for the thermal gas). As 
Fig. 3b shows, the total luminosity from the adiabatically compressed gas is much lower 
than that of the thermal gas since there are fewer high energy electrons. The slope of 
the spectrum is, however, the same for the thermal and adiabatically compressed gases 
because both are highly self-absorbed and have very similar source functions. 

The assumption of constant electron internal energy is valid at high rh where the 
cooling is efficient, as has been demonstrated by a number of previous calculations (e.g. 
Narayan & Yi, 1995b). For the lower accretion rates of interest here, where adiabatic 
compression determines the electron temperature profile, a constant electron temper- 
ature is an invalid assumption. Future calculations of synchrotron emission from low 
rh systems must self consistently calculate the evolution of the electron internal energy 
during adiabatic compression. The details of this calculation are given in Appendix B. 

Finally, we note that the slope of the synchrotron spectrum in high rh ADAFs is 
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very similar to the characteristic radio slope observed in many black hole candidates 
(Mahadevan 1997). This is, however, a direct consequence of the assumed constancy of 
9 e , which is not in general valid. Therefore we expect that low rh ADAFs will have (1) 
an unexpectedly low radio luminosity because of the absence of high energy electrons in 
the distribution function and (2) a radio slope that differs from that of high rh systems 
because the assumption of constant electron internal energy is no longer valid. 



6. Discussion & Conclusions 

In an advection dominated accretion flow, the timescale for protons and electrons to 
exchange energy by Coulomb collisions is sufficiently long compared to the accretion time 
that the two species are essentially decoupled. The protons and electrons can therefore 
have different temperatures and, if non-thermal, different distribution functions. The 
precise form of the electron and proton distribution functions is crucial for comparing 
calculated spectra from ADAFs with observations. 

Coulomb collisions among the protons are too inefficient to force the protons to be 
thermal, since the timescale for thermalization is much longer than the accretion time 
(§2). The proton distribution function is therefore determined by the viscous heating 
mechanism, which is unknown. Since the protons are marginally relativistic, Mahade- 
van et al. (1997) have shown that an ADAF produces a substantial 7-ray flux, which is 
created from the decay of neutral pions produced through proton-proton collisions. In 
particular, they have shown that the luminosity and shape of the gamma-ray spectrum 
differs dramatically if the protons have a thermal or a power-law distribution. Compari- 
son of the predicted 7-ray spectrum with observations can therefore determine the form 
of the proton distribution function and, in principle, provide a means of understanding 
the viscous heating mechanism in hot accretion flows. Mahadevan et al. (1997) have 
shown that the EGRET observations of the source 2EG J1746-2852 (Merck et al. 1996), 
which is coincident with the Galactic Center Sgr A*, are in good agreement with the 
predicted 7-ray spectrum from an ADAF, provided that the protons have a power-law 
distribution of the form E~ 2 - 75 . This is similar to the cosmic-ray proton energy dis- 
tribution, which leads to the speculation that the heating mechanism responsible for 
accelerating cosmic rays might also be at work in ADAFs. 

We have considered thermalization of electrons through both Coulomb collisions 
and self-absorbed synchrotron radiation. The latter process is significantly more impor- 
tant, since it can operate even when the electrons are effectively collisionless. We find 
that, for high accretion rates, m^4x 10~ 2 o; 2 , the electrons can efficiently exchange en- 
ergy throughout the ADAF by both Coulomb collisions and the absorption and emission 
of synchrotron photons (§3). To conclude from this that the electrons are thermal also 
requires that thermalization proceed on a shorter timescale than the heating/cooling 
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of the electrons. This ensures that the electron distribution function can relax to a 
Maxwellian more rapidly than the heating/cooling induces modifications in the distri- 
bution function. Since synchrotron radiation is the dominant cooling mechanism, and 
also leads to thermalization because it is highly self-absorbed, this criterion is trivially 
satisfied for the cooling of the electrons in ADAFs. 

Modifications of the electron distribution function by heating are perhaps a more 
viable concern. For the high mass accretion rates where thermalization appears to be 
efficient (m^4x 10~ 2 a 2 ), the electrons are primarily heated by Coulomb collisions with 
the hotter protons, provided that most of the viscously generated energy is transferred 
to the protons (§4.2). Since electron heating by the protons occurs on a timescale 
long compared to the thermalization timescale by electron-electron Coulomb collisions 
(§2.2), which is in turn long compared to the thermalization timescale by synchrotron 
self-absorption (§3.2), we conclude that heating of the electrons by the protons will 
not significantly modify the electrons from a thermal distribution.^ Therefore, for rh 
4x 10 _2 a 2 , where self-absorbed synchrotron radiation allows for efficient energy exchange 
among the electrons throughout the accretion flow, the electrons are likely to be thermal. 

We find that for lower accretion rates, rh ^ 10~ 4 a 2 , thermalization is inefficient 
throughout the ADAF and thus the electrons are not necessarily thermal. In this regime, 
the heating and cooling mechanisms determine the electron distribution function. At 
these low accretion rates, cooling is very inefficient and compressive heating of the elec- 
trons is the dominant heating mechanism (§4). The evolution of the electron distri- 
bution function is therefore determined by adiabatic compression. Using the principle 
of adiabatic invariance, we find that the electron distribution function evolves by the 
transformation p — > ap, where p is the electron momenta (relativistic or non-relativistic) 
and a is a function of the electron density. For ADAFs, in which the internal energy of 
the gas contains a contribution from the magnetic field energy, the expression for the 
compression factor, a, is more complicated than that given in §4.3. The precise form 
of a(p) does not, however, effect any of the conclusions of this paper. Nonetheless, in 
Appendix B we derive the form of a(p) under the assumption that the magnetic field 
always contributes a constant fraction of the total pressure. 

The above transformation for the electron momenta under adiabatic compression 



yields the distribution function given in equation (p8), which is non-thermal if the 
electrons have relativistic energies. In particular, the high energy tail of an adiabatically 
compressed (expanded) gas is strongly suppressed (enhanced) with respect to a thermal 



Nonthermal heating of the electrons by wave-particle interactions, reconnection, etc. could, perhaps, 
modify the electrons from a thermal distribution. We have, however, as mentioned in the introduction, 
restricted our analysis to those processes whose role in ADAFs is relatively well understood, and so have 
not considered these effects. In addition, we note that the applicability of some of these non-thermal 
heating mechanisms to the nearly collisionless plasma in ADAFs is unclear. 
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gas of the same average energy (cf. Fig. 2). We emphasize that this is not specific to 
an ADAF. Any collisionless plasma which undergoes an adiabatic change will have a 
non-thermal distribution function if the particles' momenta cross the relativistic/ non- 
relativistic boundary 

Our derivation of the distribution function of an adiabatically compressed electron 
gas is based on the assumption that the electrons are thermal at large radii. In many 
models where ADAFs have been successfully applied, the ADAF forms from a thin 
accretion disk at r ~ 10 4 (Narayan, McClintock, & Yi 1996). In this case, the electron 
distribution function is Maxwellian at large radii by virtue of the low temperature and 
large density in the thin accretion disk (cf. eq. 0). In the event that an accretion 
disk is not present, as may be the case for super massive black holes at the centers of 
elliptical galaxies which accrete the surrounding gas through Bondi accretion, the gas 
at large radii is probably still thermal. This follows from the X-ray gas profiles of the 
centers of elliptical galaxies, which indicates that the gas is thermal with a temperature 
~ 10 7 - 5 K (eg. Trinchieri et al. 1988). 

For accretion rates such that 10~ 4 ^ (m/a 2 ) ^4x 10 -2 the electrons are thermal 
for only part of the accretion flow. The initially thermal electron gas is adiabatically 
compressed down to a radius ~ rjr h (cf. eq. ||26||), with no means of thermalizing 
by synchrotron self-absorption or Coulomb collisions. For r ^ r£ h , however, electron 
thermalization by synchrotron self-absorption becomes efficient because the magnetic 
field is stronger and there is more synchrotron radiation in the plasma. 

As the accreting gas passes through the "transition" radius, r|r h , the electrons switch 
from an adiabatically compressed distribution function to a thermal one. This may 
substantially modify the predicted spectra from some ADAFs. For emission mechanisms 
in which the bulk of the emission comes from electrons with roughly the mean energy 
of the plasma, such as bremsstrahlung, the modified distribution function is unlikely to 
have a pronounced effect on the spectrum. Synchrotron radiation, on the other hand, is 
highly self-absorbed, with most of the emission coming from electrons in the high energy 
tail of the distribution function. Therefore, even when the mean electron energy is non- 
relativistic, the non-thermal distribution function due to adiabatic compression modifies 
the synchrotron spectrum substantially. In §5 we explicitly calculated the synchrotron 
spectrum from an adiabatically compressed electron gas. The synchrotron luminosity 
can be significantly smaller than that of a thermal electron gas with the same energy 
because there is a deficit of high energy particles in the tail of the distribution function 
(cf. Fig. 3). We note that, while the applications in this paper have been to ADAFs, the 
synchrotron emissivity and source function calculated in §5 apply to any adiabatically 
compressed gas. 

Previous work in calculating the spectra from ADAFs has assumed that the elec- 
trons are thermal at all radii. These systems include A0620, V404Cyg (Narayan et 
al. 1996), and NGC 4258 (Lasota et al. 1996), where the accretion rates all satisfy 
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7Ti <^ 10~ 2 a 2 . At these high accretion rates the assumption of thermal electrons is rea- 
sonably valid and the predicted spectra from these systems is unchanged by the present 
work. For systems with lower accretion rates, however, the change in the distribution 
function described here might considerably modify the predicted spectrum. 

In particular, the ADAF model of Sgr A* (Narayan et al. 1995) should be reex- 
amined since the estimated accretion rate falls below 10 _2 a 2 . Narayan et al. (1995) 
obtained good agreement with the observed spectrum of Sgr A* using m = 7 x 10 5 , 
77i ~ 3 x 10~ 3 a 2 , and an equilibrium temperature of 8 e ~ 1.5. Using these parameters 
and equation (^) gives rj; h ~ 600. For r ^ 600 the assumption of thermal electrons 
is no longer valid. Since each frequency in the radio spectrum corresponds to emission 
from a particular radius, v oc r~ 5 / 4 (Mahadevan 1997), and the maximum synchrotron 
frequency for Sgr A* is ~ 10 12 Hz (Narayan et al. 1995), the emission from r ^ 600 
corresponds to frequencies ^ 10 9 Hz. For these frequencies the model of Narayan et 
al. (1995) may be inconsistent, but since almost all of the data is for v ^ 10 9 Hz (see 
Fig. 1 in Narayan et al. 1995), the possibility of non-thermal electrons is unlikely to be 
important for this system. 

A more promising application of the present work is in modeling the 3x 10 7 M & black 
hole in M31 (Kormendy & Richstone 1995). The accretion rate in this system is esti- 
mated to be m ~ 10~ 4 (Goodman & Lee 1989), which is low enough that thermalization 
at all radii is unlikely. Using an ADAF model with a = 0.3 and T e ~ T p (r/10 4 ) 1/3 (§4.2) 
gives rf h ~ 100. For r > 100 the synchrotron spectrum would need to be calculated using 
the adiabatically compressed distribution function and would differ considerably from 
the spectrum calculated from a thermal distribution of electrons. Since lower frequencies 
in the synchrotron spectrum corresponds to emission from larger radii, the emission at 
frequencies ~ 10 10 Hz may be greatly reduced because of the new electron distribution 
function (cf. Fig. 3b). In particular, this may explain the extremely low radio flux of 
~ 10 32 ergs s _1 at v ~ 10 10 Hz (Crane et al. 1992) in M31, but detailed numerical 
calculations are required to confirm this hypothesis. 

Another potential application of this work is in determining whether elliptical galax- 
ies host dead quasars (Fabian & Canizares 1988). As suggested by Fabian & Rees (1995), 
ADAFs might allow elliptical galaxies to host ^ 10 8 M black holes and still remain 
dim. By varying the mass of the central black hole, Mahadevan (1997; see also Reynolds 
et al. 1996) has given spectra, which agree with the radio and X-ray upper limits, for 
a few of the nearby elliptical galaxies. He found that elliptical galaxies can host black 
holes up to ^ 10 9 ' 5 M & , provided that they accrete via ADAFs. Interestingly, the con- 
straint on the mass does not come from the X-ray data, but rather from the radio flux 
upper limit. Since these galaxies have estimated accretion rates (m/a 2 ) ^4x 10 _4 M 8 
(Mahadevan 1997, eqn.[56]), where M 8 = 10 8 M , the radio spectrum may be reduced 
due to the adiabatically compressed electron distribution function. This would lead to 
a further increase in the allowed mass of the central black holes, therefore completely 
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eliminating the dead quasar host problem highlighted by Fabian & Canizares (1988). 

A final application of this work is in determining accurate emission spectra from 
isolated black holes in the disk and halo of our Galaxy, which accrete from the interstellar 
medium. Ipser & Price (1982; 1977) have calculated spectra from these systems under 
the assumption that the black holes accrete by Bondi accretion and that the spectra 
is dominated by optically thin thermal synchrotron radiation. Even a small amount 
of angular momentum in the accreting gas precludes Bondi accretion, and so the gas is 
more likely to accrete onto the black hole as an ADAF. Furthermore, since these systems 
are likely to have low accretion rates (m & 10~ 5 ), the electrons are unable to thermalize 
and their distribution function is probably determined by adiabatic compression. The 
spectra from these systems will therefore be considerably different, with lower fluxes 
and luminosities peaking in different wave bands, than those obtained using a thermal 
distribution of electrons. Since the spectra from ADAFs are robust and well understood, 
detailed calculations of the spectra from isolated black holes, taking into account the 
non-thermal electron distribution function, would be quite valuable for the potential 
detection of these systems by the Sloan Digital Sky Survey. In addition, if no candidates 
are found, the ADAF spectrum will considerably modify the limits on halo black hole 
populations set by Heckler & Kolb (1996). 
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A. Asymptotic Formulae for I(x, n) and I'(x, n) 
A.l. I(x,n) 



The definition of I(x, n) is given in equation 
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This can be evaluated using the limiting forms for F(x) (Rybicki & Lightman 1979; 
eqns. [6.34a,b]), which gives 
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A. 2. J'fonl 



The function I'(x, n) is the angle averaged value of I(x, n) and is defined by 
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For i < 1 the integration is the same as done by Mahadevan et al. (1996), which gives 
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For x> 1 we use the fact that most of the emission comes from sin at p ~ 1; the integrals 
can then be evaluated by a method similar to that given in Mahadevan et al. (1996). 
We obtain 

vr 3 / 2 



I'(x,2) 
I'(x,3) 



4 

7T 3 / 2 



r~ 1/2 (l - gji/a) exp (-2x 1/2 ) ~ 1.39 £- 1/4 exp (-2 
exp (-2x 1/2 ) = 1.39 exp (-2x 1/2 ) , x > 1. 



1/2 



X > 1, 

(A6) 
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B. Compression factor in ADAFs 



As discussed is §4.3, the quantity p 3 p _1 is conserved in the adiabatic compression of 
a monatomic ideal gas, which leads to the transformation p — > ap oc p l ^p for the particle 
momenta. When the internal energy and pressure of the gas contain a contribution from 
the magnetic field, as in ADAF models, the form of the adiabatic invariant is more 
complicated. We make the assumption that the magnetic field always contributes a 
(constant) fraction /3 ac j v of the total pressure and use this to determine the expression 
for the compression factor, a. 

The gas pressure, P g , and internal energy per unit mass, U g , can be expressed in 
terms of an isotropic distribution function as 



dp 



n(p)p 4 



(p2 + 1)1/2 



(B7) 



and 



U 9 = 4n 



dpn(p)p 2 (p 2 + 1) 



1/2 



bPg/P, 



where, for simplicity, we set m = c = 1 in this Appendix. For the present purposes the 
distribution function is given by equation (|38), that is, 



n(p) 



(2tt5 2 ) 3 / 2 



exp 



— p 
~2d 



(B9) 



where a = a\/d e . We treat this distribution function and the thermodynamic quantities 
Ug and P g as functions of a and p. In this interpretation, the temperature (9 e ) which 
appears in the definition of a is a constant, namely the gas temperature when the 
adiabatic compression began. 



For a gas in an isotropically tangled magnetic field, the total pressure is 



tot 



B 2 



and the total internal energy per unit mass is 

B 2 P 



tot 



U + 



lot 



8vrp 



9 



[3 + /9adv(6 - 3)] . 



(BIO) 



(Bll) 



Adiabaticity requires that dU to t = Ptotdln pj p, which yields the following ordinary dif- 
ferential equation for the evolution of a as the gas is compressed. 



din a 
dlnp 



(3 + /3 ad v(& 



rflnPg db 

6 ,i ~ + Padv , ; Z 

dm a a In a 



1. 



(B12) 



In the non-relativistic limit, i.e., a < 1, & = 3/2, db/dlna = 0, and dlnP g /dlna 
which implies that a oc pi/( 6_3 ^dv)_ p or equipartition magnetic fields, a oc p 



2/9 K r -l/3_ 
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In the relativistic limit, a ^> 1, b — 3, db/dlna = 0, and d In P g /d In a = 1; thus 
a oc p 1 / 3 oc r^ 1 / 2 , independent of (3 rm adv We emphasize that in this context, relativistic 
vs. non-relativistic refers to the magnitude of a. It does not refer to the form of the 
distribution function, which is the same (eq. ||B9|| ) regardless of the particle energies. 

The above expressions for d(p) in the non-relativistic and relativistic limits are 
precise analogs of the expressions for T(p) for a thermal distribution given in §4.2. 
This is because the form of the distribution function (eq. |[B9|| ) is only relevant when 
calculating the evolution of a through the non-relativistic/relativistic transition (which 
must be done numerically). 



C. Departures from the Adiabatically Compressed Distribution Function 



For a nearly thermal gas, the systematic and stochastic Fokker-Planck coefficients 
given in §3 yield the timescale on which a particle's energy changes and thus the timescale 
on which the distribution becomes thermal. For non-thermal distribution functions, 
however, the distribution function may change on a timescale short compared to the 
thermalization timescale. As we now show, this implies that the adiabatically com- 
pressed distribution function given by equation Q3"8| ) is only strictly valid for smaller m 
than would be inferred from the thermalization timescales given in §3. 

We consider only self-absorbed synchrotron radiation since it is more efficient than 
Coulomb collisions in modifying the distribution function. We also assume that, even 
for the non-thermal distribution function, the Fokker-Planck coefficients are given by 
(§3.2) 

and 



(C13) 



d(A 7 y 



2dB : 



:<7>r 



(C14) 



dt ' 3 

where we have replaced the temperature of the plasma with the average 7 of the electrons, 



which is given by equation (J39f) for the adiabatically compressed distribution function. 



The timescale on which the number of particles, 71(7), at a given 7 e- folds is given 
by t n ~ l/\d\an{^)/ dt\. This, not the systematic or diffusion timescales introduced 
in §3, is the relevant timescale when considering departures from a non-thermal distri- 
bution function (rather than approaches to a thermal distribution). Evaluating t n for 
the adiabatically compressed distribution function, using the Fokker-Planck equation, 
yields (cf. eqs. 0, |3§) 



t , 



1 



7 



1047 2 



+ 



+ 



64 7 4 



(7) 37r( 7 ) 2 7r( 7 ) 3 3vr 2 (7) 



(C15) 
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where we have assumed that the electrons are relativistic. For large 7, the timescale 
on which the distribution function changes is ~ {l/(l)) A times faster than the usual 
diffusion time, t s diS , and ~ (7/ (7)) 3 times faster than the usual systematic timescale, tl ys 
(eq. PBfl). The physical reason for this is that diffusion can be much more efficient than 
is indicated by the stochastic Fokker-Planck coefficient if the non-thermal distribution 
function varies strongly with energy. 

For accretion times less than the minimum of the diffusion and systematic timescales 
(t a < min[t diff , tg ys ]) it is incorrect to assume a thermal distribution function, but it is 
only if t a < t n that it is strictly valid to assume that adiabatic compression yields the 
distribution function given by equation (|38|). This occurs provided that 

. < 7.9 x IP' 5 .7 1047 2 8 7 3 64 7 4 -1 

— ; ra 8 — 4—— H oH r (C16) 

(7) V 0-5 / (7) 3tt( 7 ) 2 tt( 7 ) 3 3tt 2 ( 7 ) 4 1 ^ 

For 7 pa (7) there are no strong gradients in the distribution function and so the dis- 
tribution function changes only on the energy diffusion time, i.e., t n ~ t s diS . However, 
for 7 ph 3(7), which is the Lorentz factor of the electrons responsible for most of the 
synchrotron emission, equation ( |C16| ) shows that the electrons have an adiabatically 
compressed distribution function only for accretion rates ^ 10~ 6 a: 2 r. This is a substan- 
tially more stringent requirement than the condition that thermalization be inefficient. 
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Figure Captions. 



Figure 1: For accretion rates above a given curve in the figure, an electron with Lorentz 
factor 7 is able to thermalize with a background plasma of temperature 9 e . The dashed 
lines represent thermalization through Coulomb collisions, while the solid lines repre- 
sent thermalization through synchrotron self-absorption. The curves are determined by 
whether the dominant thermalizing process is systematic or diffusive acceleration (cf. 
eqs. ||15|| , ||25|| ). While thermalization by Coulomb collisions is, for a fixed temperature, 



independent of the radius r (cf. eq.|T5[), the accretion rates required for thermalization 



by synchrotron self-absorption increase linearly with r (cf. eq. f25|). 

Figure 2: a) A comparison of the momentum distribution function for a thermal gas 
of temperature 9 e ~ 1 with that of an adiabatically compressed collisionless gas of the 
same energy. The adiabatically compressed gas is initially thermal and non-relativistic, 
but is then compressed to relativistic energies, resulting in a non-thermal distribution 
function, b) Analogous to a), but for the adiabatic expansion of a gas from relativistic 
to non-relativistic energies. 

Figure 3: a) A comparison of the synchrotron radiation from a sphere of radius ~ 7x 10 11 
cm for an adiabatically compressed (solid line) and thermal (dashed line) gas with the 
same average energy. The electron number density and magnetic field strength are 
constant throughout the sphere (n e ~ 10 10 cm -3 and B ~ 200 Gauss). At low frequen- 
cies the radiation is highly self-absorbed, and the observed intensity is given by the 
source function, not the emissivity (dotted line) of the respective gases, b) The total 
synchrotron spectrum from an ADAF (m = 2 x 10 6 , rh = 10~ 5 ) for adiabatically com- 
pressed electrons (solid line) and thermal electrons (dashed line) of the same average 
energy, which is assumed to be constant throughout the ADAF. This spectrum is ob- 
tained by summing the individual spectrum from each radius in the accretion flow, and 
thus accounts for the variation in electron number density and magnetic field strength 
with radius. 
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Figure 2 
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Figure 3 
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